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ABSTRACT 

It has been recently shown (Munyaneza, Tsiklauri and Viollier 1998) that the analysis of the 
orbits of the fast moving stars close to Sgr A* (Eckart and Genzel 1997) provides a valuable 
tool to probe the gravitational potential near the Galactic center. As an example, we present 
here the results on a calculation of possible orbits of the star SO-1 in both, the black hole and 
degenerate neutrino ball scenarios of the central mass, based on the recent measurements of 
stellar proper motions at the Galactic center by Ghez et al. 1998. Taking into account the error 
bars of their analysis, it is shown that within a few years time, the orbit of SO-1 may indeed 
reveal the nature of the supermassive compact dark object at the Galactic center. 



Subject headings: black hole physics — celestial mechanics, stellar dynamics — dark matter 
elementary particles — Galaxy: center 



1. Introduction 

The determination of the mass distribution near the center of our Galaxy and the question, whether 
it harbours a supermassive black hole (BH) or not, have been long-standing issues (Oort 1977; Genzel and 
Townes 1987; Genzel et al. 1994 and Ho 1998 for a recent review). Various techniques have been used to 
find the mass of this supermassive compact dark object which is usually identified with the radio source 
Sagittarius A* (Sgr A*) at or near the Galactic center. The most detailed information to date comes 
from the statistical analysis of the dynamics of stars moving in the gravitational field of the central mass 
distribution (Sellgren et al. 1987; Rieke and Rieke 1988; McGinn et al. 1989; Sellgren et al. 1990; Lindqvist 
et al. 1992; Haller et al. 1996; Eckart and Genzel 1996; Genzel et al. 1996; Eckart and Genzel 1997; Genzel 
et al. 1997; Ghez et al. 1998). Genzel et al. 1997 have established that the central dark object has a mass 
of (2.61 ± 0.76) x 1O 6 M , concentrated within a radius of 0.016 pc and located very close to Sgr A*. In the 
most recent observations, Ghez at al. 1998 confirm a mass of (2.6 ± 0.2) x 10 6 M Q , enclosed within a radius 
of 0.015 pc. In the latter observations, the accuracy of the velocity measurements in the central arcsec 2 has 
been improved considerably, and thus the error bar on the central mass has been reduced by about a factor 
of 4. In both data sets, the presence of a supermassive compact dark object is revealed by the fact that 
several stars are moving within a projected distance of less than 0.01 pc from the central radio source Sgr 
A* at projected velocities in excess of 1000 km/s. 

For completeness, we mention here that the mass distribution at the Galactic center could also be 
studied through the motion of gas clouds and streamers (Lacy et al. 1980; Genzel & Townes 1987; Lacy et 
al. 1991). However, gas flows may be easily perturbed by non-gravitational forces such as shocks, radiation 
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pressure, winds, magnetic fields, etc., and hence this probe is considered to be less reliable for determining 
the mass of the compact dark object at the Galactic center. 

The non-thermal spectrum of Sgr A* (Serabyn et al. 1997), that has been shown to originate from 
a very compact source (Rogers et al. 1994; Genzel et al. 1997; Ghez et al. 1998), and the low proper 
motion of Sgr A* (Backer 1996) have led many (e.g. Lynden-Bell and Rees 1971) to suggest that Sgr A* 
may be a supermassive BH of mass ~ 2.6 x 1O 6 M0. Supermassive BHs have also been inferred for several 
other galaxies such as M87 (Ford et al. 1994; Harms ct al. 1994; Macchetto et al. 1997) and NGC 4258 
(Greenhill et al. 1995; Myoshi et al. 1995). Taking this suggestion seriously, one is immediately faced 
with fundamental issues such as the prevalence of supermassive BHs in the nuclei of normal galaxies and 
the nature of the accretion mechanism that makes Sgr A* so much fainter than typical active galactic 
nuclei (Melia 1994; Narayan et al. 1995). However, as the best current observations probe the gravitational 
potential at radii 4 x 10 4 larger than the Schwarzschild radius of a BH of mass 2.6 x 10 6 M Q (Ghez et al. 
1998), it is perhaps prudent not to focus too much on the BH scenario, without having explored alternative 
scenarios for the supermassive compact dark object. 

One alternative to the BH scenario is a very compact stellar cluster (Haller et al. 1996, Sanders, 1992). 
However, based on the evaporation and collision time stability criteria, it is doubtful that such clusters 
could have survived up to the present time (see Moffat 1997 for an alternative point of view). Indeed, in the 
case of our Galaxy and NGC 4258, Maoz (1995,1998) has found that even the lower limits to the half-mass 
densities of such compact clusters (1 x lO 12 M pc~ 3 for NGC 4258 and 6 x 10 n M Q pc~ 3 for our Galaxy) 
are too large that they could be due to stable clusters of stellar or substellar remnants. The estimated 
maximal lifetimes for such dense clusters are about 10 s years for our Galaxy and a few 10 8 years for the 
NGC 4258, i.e. much shorter than the age of the Universe. This seems to rule out the existence of dense 
clusters at the centers of the above mentioned galaxies, unless we are prepared to believe that we happen 
to live in a privileged epoch of the lifetime of the Universe. Note, however, that for other galaxies, such 
as M31, M32, M87, NGC 3115, NGC 3377, NGC 4261, NGC 4342, NGC 4486B and NGC 4594, maximal 
lifetimes of dense stellar clusters are in excess of 10 11 years. Moreover, it should be acknowledged that the 
uncertainties in the understanding of the core collapse process of such dense clusters still leave some room 
for speculation about a possible interpretation of the supermassive compact dark objects at the centers of 
galaxies (including both, our Galaxy and NGC 4258) in terms of e.g. core-collapsed clusters (Maoz 1998). 
But, apart from a cluster of very low mass BH's that is free of stability problems, the most attractive 
alternative to a dense stellar cluster is a cluster of elementary particles. 

In fact, in the recent past, an alternative model for the supermassive compact dark objects in galactic 
centers has been developed (Viollier et al. 1992, 1993; Viollier 1994; Tsiklauri and Viollicr 1996, 1998a,b, 
1999; Bilic et al. 1998; Bilic et al. 1999). The cornerstone of this model is that the dark matter at 
the center of galaxies is made of nonbaryonic matter in the form of massive neutrinos that interact 
gravitationally forming supermassive neutrino balls in which the degeneracy pressure of the neutrinos 
balances their self-gravity. Such neutrino balls could have been formed in the early Universe during a 
first-order gravitational phase transition (Bilic and Viollier 1997, 1998, 1999a, b). In fact, it has been recently 
shown that the dark matter concentration observed through stellar motion at the Galactic center (Eckart 
& Genzel 1997; Genzel et al. 1996) is consistent with a supermassive object of 2.5 x 10 6 solar masses 
made of self-gravitating, degenerate heavy neutrino matter (Tsiklauri & Viollier 1998a). Moreover, it has 
been shown that an acceptable fit to the infrared and radio spectrum above 20 GHz , which is presumably 
emitted by the compact dark object, can be reproduced in the framework of standard accretion disk theory 
(Tsiklauri & Viollier 1999; Bilic et al. 1998), in terms of a baryonic disk immersed in the shallow potential 
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of the degenerate neutrino ball of 2.5 x 10 6 solar masses. 

The purpose of this paper is to compare the predictions of these two models for the supermassive 
compact dark object at the center our Galaxy, i.e. (i) the black hole scenario and (ii) the degenerate 
neutrino ball scenario as an example of an extended object. Both these models are not in contradiction 
with the technologically challenging proper motions observations and their statistical interpretation (Gcnzel 
et al. 1997; Ghez et al. 1998) . It is therefore desirable to have an additional independent dynamical 
test, in order to distinguish between these two possible scenarios describing the compact dark object at 
the center of our Galaxy. In the recent past, mainly statistical arguments involving many stars have been 
used to determine the gravitational potential at the Galactic center. However, in this paper, we would like 
to demonstrate that it is also possible to draw definite conclusions from the motion of individual stars, in 
particular, in the immediate vicinity of the Galactic center, where statistical arguments cannot be easily 
applied due to the low density of stars. To this end, we have recently calculated the orbits (Munyaneza , 
Tsiklauri and Viollier 1998) of the fastest moving infrared source SI using the Genzel et al. 1997 data for 
a supermassive BH or a neutrino ball mass of 2.61 x 10 6 solar masses. We have shown that tracking the 
orbits of SI offers a good opportunity to distinguish in a few years time between the two scenarios for the 
supermassive compact dark object. Here we perform a full analysis of the orbits of the same star SO-1 based 
on the most recent Ghez et al. 1998 data, including all the error bars of the measurements. A distance to 
the Galactic center of 8 kpc has been assumed throughout this paper. 

This paper is organized as follows: In section 2, we present the equations that describe degenerate 
neutrino balls and we establish some constraints on the neutrino mass based on Ghez et al. 1998 data. In 
section 3, we study the dynamics of SO-1 and conclude with the discussion in section 4. 



2. The compact dark object as a neutrino ball 

Dark matter at the Galactic center can be described by the gravitational potential $(r) of the neutrinos 
and antincutrinos that satisfies Poisson's equation 

A$ = 4^Gp„ (1) 

where G is Newton's gravitational constant and p v is the mass density of the neutrinos and antincutrinos. 
Neutrino matter will interact gravitationally to form supermassive neutrino balls in which self-gravity of 
the neutrinos is being balanced by their degeneracy pressure P v (r) according to the equation of hydrostatic 
equilibrium 

dP v d<P 

-7- = -Pv~T ■ ( 2 ) 

dr dr 

In order to solve equation (|l|) , one needs a relation between the pressure P v and the density p v . To this 
end we choose the polytropic equation of state of degenerate neutrino matter, i.e. 

Pv = Kpl'\ (3) 

where the polytropic constant K is given by (Viollier, 1994) 

\yv J bmj 

Here, m„ denotes the neutrino mass, g v is the spin degeneracy factor of the neutrinos and antincutrinos, i.e. 
g u = 2 for Majorana and g u — 4 for Dirac neutrinos and antincutrinos. We now introduce the dimensionless 
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potential and radial variable, v and x , by 



GM Q ( , 

$(r) = — w'(x ) - — , (5) 



r = a u x, (6) 

where xq is the dimensionless radius of the neutrino ball, and the scale factor a v which plays here the role 
of a length unit is given by 

a „= 2 , 13761yrx (^) !, V/, (7) 

Assuming spherical symmetry, we finally arrive at the non-linear Lane-Emden equation 

d 2 v v 3 / 2 



dx 2 x 1 / 2 



(8) 



with polytropic index 3/2. The boundary conditions are chosen in such a way that v vanishes at the 
boundary xo of the neutrino ball. The mass M B of a (pointlike) baryonic star at the center of the neutrino 
ball is fixed by v(0) = M B /M Q . The case M B — corresponds to a pure neutrino ball without a pointlikc 
source at the center. The mass enclosed within a radius r in a pure neutrino ball can be written in terms of 
v{x) and its derivative v'(x) as 

M(r) = / ^p v r 2 dr = -M© (v'(x)x - v{x)) . (9) 



In order to describe the compact dark object at the Galactic center as a neutrino ball and constrain 
its physical parameters appropriately, it is worthwhile to use the most recent observational data by Ghez et 
al. 1998, who established that the mass enclosed within 0.015 pc at the Galactic center is (2.6 ± 0.2) x 10 6 
solar masses. Following the analysis of Tsiklauri & Viollier 1998a, we choose the minimal neutrino mass 
m v to reproduce the observed matter distribution, as can be seen from Fig. 1, where we have added 
the Ghez et al. 1998 and Genzel ct al. 1997 data points with error bars. In Fig. 1 we include only the 
neutrino ball contribution to the enclosed mass, as the stellar cluster contribution is negligible by orders 
of magnitude at these radii. For a M = 2.4 x 10 6 M Q neutrino ball, the constraints on the neutrino mass 
are vn v > 17.50 keV/c 2 for g v = 2 and m v > 14.72 keV/c 2 for g v — 4, and the radius of the neutrino 
ball is R < 1.50 x 10~ 2 pc. Using the value of M = 2.6 x 10 6 M Q , the bounds on the neutrino mass are 
m„ > 15.92 keV/c 2 for g v — 2 or m v > 13.39 keV/c 2 for g v = 4 and the radius of the neutrino ball turns 
out to be R < 1.88 x 10~ 2 pc . Finally, for a M = 2.8 x 10 6 M Q neutrino ball, the range of neutrino mass 
is m v > 15.31 keV/c 2 for g v = 2 and vn v > 12.87 keV/c 2 for g v — 4 and the corresponding neutrino ball 
radius R < 2.04 x 10~ 2 pc. 



3. Dynamics of SO-1 

We investigate the motion of S0-1 that is the star closest to the Galactic center, and at the same 
time, also the fastest of the 15 stars in the central arcsec 2 around Sgr A*. We study the motion of S0-1 
in the gravitational potential near Sgr A*, assuming that the central object is cither a BH of mass M 
or a spatially extended object represented by a neutrino ball of mass M, that consists of self-gravitating 
degenerate heavy neutrino matter. The BH or neutrino ball mass M will be taken to be 2.4, 2.6 and 2.8 
x 10 6 solar masses which corresponds to the range allowed by the Ghez et al. 1998 data. We use Newtonian 
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dynamics, as the problem is essentially nonrelativistic, because the mass of the neutrino ball is much less 
than the Oppenheimer-Volkoff limit corresponding to this particular neutrino mass (Bilic, Munyaneza & 
Viollier 1999). Consequently, we can write Newton's equations of motion as 

GM(r) ^ 



(x 2 + y 2 + z 2 ) 3 / 2 
GM(r) 



y 2 + z 2 fl 2 ' 



(11) 



■ z= z m 

(x 2 + y 2 + z 2 fl 2 ' 1 ' 



where x, y, z denote the components of the radius vector of the star SO-1 and r = \J x 2 + y 2 + z 2 , Sgr A* 
being the origin of the coordinate system. We thus assume that the center of the neutrino ball and the BH 
is at the position of Sgr A*. The dot denotes of course the derivative with respect to time. In the case of 
a BH, M (r) = M is independent of r, while in the neutrino ball scenario, M(r) is given by Eq. (||) and it 
reaches M(R) = M at the radius of the neutrino ball R. The initial positions and velocities for this system 
of equations are taken to be those of SO-1 in 1995.4, when the coordinates of SO-1 were RA = —0.107" 
and DEC = 0.039". The x and y components of the projected velocity are v x = 470 ± 130 km/s and 
v y = —1330 ± 140 km/s (Chez et al. 1998), respectively. Here x is opposite to the RA direction and y is in 
the DEC direction. 

In Fig. 2 we plot two typical orbits of S0-1 corresponding to a BH and neutrino ball mass of 
M = 2.6 x 10 6 M©. The input values for v x and v y are 470 km/s and -1330 km/s, respectively. The 
z-coordinate of the star S0-1 is assumed to be zero and the velocity component in the line-of-sight of the 
star S0-1, v Z) has also been set equal to zero in this graph. The filled square labels denote the time in years 
from 1990 till 2015. In the case of a BH, the orbit of S0-1 is an ellipse, with Sgr A* being located in one 
focus (denoted by the star in the figure). The period of S0-1 is 12.7 years and the minimal and maximal 
distances from Sgr A* are 1.49 and 7.18 light days, respectively. In the case of a neutrino ball, the orbit will 
be bound but not closed, with minimal and maximal distances from Sgr A* of 3.98 and 42.07 light days, 
respectively. It can be seen from Fig. 2 that, in the case of a neutrino ball, S0-1 is deflected much less than 
for a BH, as the gravitational force at a given distance from Sgr A* is determined by the mass enclosed 
within this distance. Using Eq. (||) we can estimate the mass enclosed within a radius corresponding to 
the projected distance of S0-1 from Sgr A* (4.41 x 10~ 3 pc ) to be ~ 1.8 x 10 5 M©. Thus, in the case of 
a neutrino ball, the force acting on S0-1 is about 14 times less than in the case of a BH. This graph can 
serve to establish, whether Sgr A* is a BH or an extended object, due to the fact that the positions of S0-1 
will differ as time goes on in the two scenarios. However, this conclusion is perhaps too optimistic, as we 
have not yet considered (i) the uncertainties in v x and v y , (ii) the error bars in the total mass of the BH or 
neutrino ball, (iii) the complete lack of information on z and v z . 

As a next step, we investigate the dependence of the orbits on the uncertainties in the velocity 
components. The results of this calculation are presented in Fig. 3 where we have set z — v z — 0. In the 
case of a BH, the orbits of S0-1 are ellipses , while the other 5 thick lines are bound orbits of S0-1 for the 
neutrino ball scenario. The spread of the orbits induced by the error bars in v x and v y is small compared 
to that of the recent analysis based on the Genzel et al. 1997 data (Munyaneza, Tsiklauri & Viollier 1998). 
The time labels, represented by filled squares on the orbits, are placed in intervals of 5 years: starting from 
1995.4 up to 2005 in the case of a BH, and up to 2015 in the case of a neutrino ball. The periods of S0-1 
for different orbits vary between 10 and 17 years for the BH scenario. We thus see that the error bars in v x 
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and v y do not alter the predictions of Fig. 1 in substance. We now would like to study, how the orbits are 
changed if we let the mass of the neutrino ball or the BH vary within the estimated error bars (Ghez et al. 
1998). 

In Fig. 4 and 5, we present the results of our calculations, for both scenarios, with central masses of 
M = 2.4 x 10 6 M Q and M = 2.8x 10 6 M Q , respectively. The neutrino masses consistent with the Ghez et al. 
(1998) data are are m v > 17.50 keV/c 2 for a M = 2.4 x 1O 6 M neutrino ball and m v > 15.31 keV/c 2 for a 
M = 2.8 x 1O 6 M neutrino ball. The filled squares represent the time labels spaced by 5 year intervals as in 
Fig. 3. In the BH scenario, the periods of S0-1 with M = 2.4 x 10 6 M Q vary between 11 years and 20 years, 
while in the case of a M — 2.8 x 10 6 M Q , the periods vary between 9.5 and 15 years. Comparing the orbits 
of S0-1 in Fig. 4 and 5 with those in Fig. 3, we conclude that the errors bars in the total mass of the BH or 
neutrino ball make no qualitative difference for the motion of S0-1. In both scenarios of the supermassive 
compact dark object, all the orbits considered for three different values of the BH or neutrino ball mass 
are bound for z — v z — 0, as can be seen from Fig. 6 and 7, where the escape and circular velocities are 
plotted as functions of the distance from Sgr A*. In these graphs, we have also included the Ghez et al. 
(1998) data with error bars, for the 15 stars in the central arsec 2 , assuming that the velocity component 
and distance from Sgr A* in the line-of-sight are both zero, i.e. v z = and z = 0. Thus, the data points are 
lower bounds on the true circular or escape velocity and radius, and the real values lie in the quarter-plane 
to the right-and-up of the measured data point. For instance, the innermost data point describing the star 
S0-1 is in both scenarios, consistent with a bound orbit if \z\ and \v z \ are not too large, as can be seen from 
the escape velocity in Fig. 6. However, S0-1 cannot be interpreted as a virialized star in the neutrino ball 
scenario, as is evident from the plot of the circular velocity in Fig 7; it thus would have to be an intruder 
star. If the projected velocity of a star at a given projected distance from Sgr A* is larger than the escape 
velocity at the same distance (assuming z = 0), the neutrino ball scenario is virtually ruled out, since the 
kinetic energy of the star would have to be very large at infinity. 

We now turn to the investigation of the dependence of the orbits on the z coordinate and z component 
of the velocity of the star S0-1. The two quantities, z and v z , are the major source of uncertainty in 
determining the exact orbit of the star S0-1. However, this shortcoming will not substantially affect the 
predictive power of our model, as we will see below. In Fig. 8 we show the results of a calculation of the 
dependence of the orbit on z for a M = 2.6 x 10 6 M Q neutrino ball or BH. The input values for v x and 
v v are fixed at 470 km/s and -1330 km/s, respectively, and v z is assumed to be zero. The z-coordinate is 
varied from zero up to the radius of the neutrino ball, i.e. the distance from Sgr A* beyond which there 
is obviously no difference between the BH and the neutrino ball scenarios. In this case, the radius of the 
neutrino ball 1.88 x 10~ 2 pc or 0.485 ". The top panel represents the orbits in the case of a BH, for different 
values of z, while the lower panel describes the dependence of the orbit on z in the neutrino ball scenario. 
We conclude from this plot that, increasing \z\ has the effect of shifting the orbits towards the lower right 
corner of the graph. This is, obviously, due to the fact that increasing \z\ means going further away from 
the scattering center, thus yielding less deflection of the orbit. Moreover, in the neutrino ball scenario, the 
dependence on z is relatively insignificant, as long as \z\ is smaller than the radius of the neutrino ball. 
This is in accordance with the fact that for small distances from the center, the potential of a neutrino ball 
can be approximated by a harmonic oscillator-type potential, where the Newtonian equations of motion 
decouple in Cartesian coordinates. The dependence of the orbits of S0-1 on v z has a similar effect as in the 
previous graph. Here, we have fixed z to zero and v z has been varied as an input parameter. Increasing \v z \ 
yields a greater velocity of the star and, obviously, a fast moving star will be deflected less than a star with 
smaller \v z \. The results of this calculation are summarized in Fig. 9. 



- 7- 



4. Conclusion and discussion 

We have demonstrated that the orbits of SO-1 differ substantially for the BH and neutrino ball scenarios 
of the Galactic center, especially with the new Ghez et al. (1998) data. We have shown that using these 
data, the error bars in velocities of SO-1 and mass of the central object do not change the pattern of the 
orbits of SO-1. In the case of a BH, the orbit of SO-1 is much more curved than in the neutrino ball scenario, 
as long as \z\ is smaller than the radius of the neutrino ball. Increasing \z\ and \v z \ shifts the orbits to the 
lower right corner of the graph, and this gives us a key to establish the allowed regions of SO-1 depending on 
whether it is a BH or a neutrino ball, irrespective of the values of the parameters z and v z . In Fig. 10 three 
orbits are plotted: the upper-leftmost orbit of SO-1 corresponding to the neutrino ball scenario (actually, 
line 9 of Fig. 4) and two orbits in a BH scenario with the smallest minimal and maximal distances from Sgr 
A* (ellipses 2 and 4 from Fig. 5). This figure serves as a test to distinguish the supermassive BH scenario 
from the neutrino ball model of the Galactic center. It is clear that, as the observations proceed within the 
next year, one might be able to tell the difference between the two models of the supermassive compact 
dark object at the center of our Galaxy. 

If the star is found in the region F inside the ellipses, this will rule out both the BH and the neutrino 
ball scenario of Sgr A*, as seen in Fig. 10. We can estimate the minimal distance of approach to Sgr A* to 
be 0.909 light days. If the orbit of SO-1 ends up in the upper-left zone of the thick line, this will clearly rule 
out the neutrino ball scenario for the chosen lower limit of the neutrino mass. However, if SO-1 is found 
in the lower right corner of the same line (i.e. below the thick line), then the supermassive object can be 
interpreted as either a neutrino ball or a BH with a large z or v z parameter. One can of course repeat this 
analysis for several stars in the central arcsec 2 and use some statistical arguments: should there be no stars 
in the black hole zone, and many stars found in the zone for black holes and neutrino balls, the black hole 
interpretation of the supermassive compact dark object at the Galactic center would become less attractive, 
as some of the stars should be moving close to the plane perpendicular to the line-of-sight, i.e. they should 
have small \v z \ and \z\. 

The neutrino masses used for the neutrino ball are lower limits. We note that increasing the neutrino 
mass will make the radius smaller (the neutrino ball radius scales as oc m^ 8 ^ 3 ) and, when it reaches the mass 
corresponding to the Oppenheimer-Volkoff limit, there will be little difference between the two scenarios. 
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Figure captions: 

Fig. 1: The mass enclosed within a distance of the center of a neutrino ball of 2.4, 2.6, 2.61 and 2.8 
millions solar masses. Based on the Ghcz ct al. 1998 data, the bounds on the neutrino mass are m v > 17.50 
keV/c 2 for g v = 2 or m v > 14.72 keV/c 2 for g v = 4 and M = 2.4 x 10 6 M Q . For M = 2.6 x 10 6 M Q the 
bounds on the neutrino mass are m v > 15.92 keV/c 2 for g v = 2 and m v > 13.39 keV/c 2 for g v = 4. Finally, 
a neutrino mass range of m„ > 15.31 keV/c 2 for g v = 2 or m v > 12.87 keV/c 2 for g u = 4 is consistent with 
a supermassive object of M = 2.8 x 10 6 M Q . The Ghez et al. 1998 and Genzel et al. 1997 data points with 
error bars are also shown in this graph. 

Fig. 2: Projected orbits of the star S0-1 for BH and neutrino ball scenarios with M = 2.6 x 10 6 M Q 
and v z = z — 0. The velocity components of S0-1 are taken to be v x = 470 km/s and v y = —1330 km/s. 
The filled squares denote the time labels. The period of S0-1 in a the BH scenario is 12.7 years and the 
minimal and maximal distances from Sgr A* are 1.49 and 7.18 light days. The orbit of S0-1 in the neutrino 
ball scenario is bound with minimal and maximal distances from Sgr A* of 3.98 and 42.07 light-days, 
respectively. 

Fig. 3: Projected orbits of the star S0-1 in the case of a BH or a neutrino ball of M = 2.6 x 1O 6 M 
taking into account the error bars in the velocity components. The labels for the different orbits are: 1: 
v x = 470 km/s and v y = —1330 km/s (median values), 2: v x = 340 km/s and v y = —1190 km/s, 3: v x = 340 
km/s and v y = —1470 km/s, 4: v x = 600 km/s and v y = —1190 km/s, 5: v x = 600 km/s and v y = —1470 
km/s. The periods of S0-1 for different orbits in the BH scenario vary between 10 and 17 years. The 
thick lines 6 to 10 correspond to the orbits in the neutrino ball scenario with the following description: 6: 
v x = 470 km/s and v y = —1330 km/s (median values), 7: v x = 340 km/s and v y — —1190 km/s, 8: v x = 340 
km/s and v y = —1470 km/s, 9: v x = 600 km/s and v y = —1190 km/s, 10: v x = 600 km/s and v y — —1470 
km/s. All the orbits in both scenarios are bound. The time labels (filled squares) on the orbits are placed 
in intervals of 5 years, up to the year 2005 in the case of a BH and up to 2015 in the case of a neutrino ball. 

Fig. 4: Projected orbits of the star S0-1 in the case of a BH or neutrino ball with M = 2.4x 1O 6 M . In 
this graph we explore how the orbits are affected by the uncertainty in the mass of the BH or neutrino ball. 
The orbits are calculated for z = v z = . The description of the orbits are the same as Fig. 3. The periods 
of S0-1 in the BH scenario vary between 11 and 20 years and all the orbits are bound in both scenarios. 

Fig. 5: Projected orbits of the star S0-1 for M — 2.8 x 10 6 M Q . All the orbits are bound and calculated 
for different values of the velocity components as in Fig. 3. The periods of S0-1 vary between 9.5 years and 
14.7 years in the BH scenario. 

Fig. 6: The escape velocity as a function of the distance from Sgr A* for BH and neutrino ball scenarios. 
The value of the mass of the central object is varied as indicated on the graph. The data points with error 
bars of 15 stars in the central arcsec 2 are taken from Ghcz ct al. 1998 assuming that the projected velocity 
and distance from Sgr A* are equal to the true velocity and distance, respectively, i.e. z = and v z = 0. 
This graph shows that S0-1 is bound in both scenarios for different values of the mass of the central object. 

Fig. 7: The circular velocity as a function of the distance from Sgr A* for BH and neutrino ball 
scenarios. The mass of the central object is varied as indicated on the graph. The data points with error 
bars of 15 stars in the central arcsec 2 are taken from Ghcz et al. 1998 assuming that the projected velocity 
and distance from Sgr A* are equal to the true velocity and distance, respectively, i.e. z = and v z = 0. 
This graph shows that the orbits of S0-1 arc almost circular in the case of the BH scenario (see text for the 
discussion concerning the innermost data point). 
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Fig. 8: Projected orbits of the star SO-1 in the case of a supcrmassive BH (top panel) and in the case of 
a neutrino ball (lower panel) with M = 2.6 x 10 6 M©. In this graph we explore how the orbits are affected 
by the uncertainty in the z-parametcr. The labels for the orbits are given in the graph. Note, that for 
z = 0.4849", which corresponds to the radius of the neutrino ball for the assumed distance to the Galactic 
center, the orbits for a BH and neutrino ball are identical, as it should be. In this graph v x = 470 km/s, 
v y = —1330 km/s and v z = 0. 

Fig. 9: Projected orbits of the star SO-1 in the case of a supermassive BH (top panel) and in the case of 
a neutrino ball (lower panel) of M = 2.6 x 10 6 M Q . In this graph we explore how the orbits are affected by 
the uncertainty in v z . The labels for the orbits are given in the graph. Here, v x = 470 km/s, v y = —1330 
km/s and z = 0. 

Fig. 10: Prediction regions for the supcrmassive central object. This graph combines line 9 from Fig. 4 
and lines 2 and 4 from Fig. 5. If the star S0-1 will be found inside the ellipses (region F), this will rule 
out both the BH and the neutrino ball models . If the star S0-1 will eventually be found in the upper- left 
zone of the graph, i.e. up and left of the thick orbit, this will rule out the neutrino ball interpretation for 
the chosen neutrino mass. Finally, if S0-1 will be found to the right and below the thick line, then the 
supermassive central object should be interpreted either as a BH with large z or as a neutrino ball. 
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